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THEORY  OF  PLASTIC  STRUCTURES 

1.  Introduction  When  a tension  specimen  is  loaded 
monotonically , it  will  generally  have  an  initial  linear  range, 
followed  by  a nonlinear  one.  Figure  1,  for  example,  shows 
experimental  results  for  mild  steel  [1]  and  245-7  aluminum  [2] . 
Figure  lb  also  shows  two  fitted  curves  for  the  aluminum  data. 
Although  the  solid  curve  is  certainly  the  better  fit,  the 
piecewise  linear  dashed  curve  has  obvious  theoretical  advantages. 
For  simplicity  of  exposition  we  shall  assume  that  it  adequately 
represents  the  "real"  material.  Points  on  the  first,  steeper, 
portion  of  the  curve  are  called  "elastic",  and  points  on  the  other 
portion  are  referred  to  as  "plastic".  A structure  is  called 
elastic  when  all  of  its  parts  correspond  to  the  elastic  portion 
of  the  curve.  A theory  of  plastic  structures  must  be  used  when 
any  part  of  the  structure  is  on  the  plastic  portion. 

It  is  convenient  to  define  several  idealized  models  of 
material  behavior  as  indicated  in  Fig.  2.  The  salient  features 
of  these  models  can  be  illustrated  by  the  simple  3-bar  truss  in 
Fig.  3.  The  results  of  loading  the  truss  with  a slowly,  mono- 
tonically increasing  load  are  shown  in  Fig.  4.  The  behavior  of 
the  "real"  truss,  curve  (a) , may  be  characterized  as  follows. 

For  P < P^,  the  elastic  limit,  the  truss  is  in  the  elastic  range— 
all  three  bars  are  elastic.  Between  P^  and  Pg,  the  yield-point 
load,  the  central  bar  is  plastic,  but  the  two  diagonal  bars  are 
still  elastic.  Since  these  bars  by  themselves  would  constitute 
a load-carrying  structure,  the  response  of  the  3-bar  truss  is 
still  of  an  elastic  order  of  magnitude  and  the  truss  is  in  the 


range  of  contained  plastic  deformation.  Along  BC,  all  three 
bars  are  plastic  and  a small  increment  of  load  produces  a much 
larger  increment  of  deformation,  so  that  we  refer  to  the  range  of 
extensive  plastic  deformation. 

We  consider  now,  the  various  idealizations  of  Fig.  3.  The 
fully  elastic  model  (b)  is  exact  in  the  elastic  range,  an 
increasingly  poor  approximation  in  the  range  of  contained  plastic 
deformation,  and  an  absurd  approximation  for  P > P^.  The  rigid/ 
strain-hardening  model  (c)  does  not  distinguish  between  the 
elastic  and  contained  ranges,  but  is  a good  approximation  to  the 
extensive  plastic  deformation  range.  The  elastic/perfectly- 
plastic  model  (d)  is  exact  in  the  elastic  range  and  a very  good 
approximation  in  the  range  of  contained  plastic  deformation. 

It  is  also  essentially  exact  for  moderate  strains  for  the  mild 
steel  of  Fig.  la.  However,  according  to  this  model,  the  truss 
cannot  support  any  load  greater  than  the  yield-point  load 


Pq  = (1  + /2)A  Oq 


(1.1) 


so  that  it  is  more  appropriate  to  refer  to  the  horizontal  segment 
of  the  curve  as  unrestricted  plastic  flow.  Any  load  greater  than 


Pq  could  not  be  equilibrated  and  hence  would  produce  accelerated 
plastic  flow.  Finally,  the  simple  rigid/perfectly-plastic  model 
(a)  does  not  give  a good  approximation  in  any  range,  but  does 
predict  the  yield-point  lead  (1.1) . 

Indeed,  all  three  of  the  idealized  plastic  models  have  exactly 
the  scune  yield-point  load  Pq,  and  this  value  is  a very  good 
approximation  to  Pg  for  the  real  material.  Therefore,  the  simple 
rigid/perfectly-plastic  model  can  be  used  to  predict  this 


significant  property  of  the  structure, 


Based  upon  this  simple  illustration,  we  may  propose  the 
following  approach  to  structural  analysis.  First,  use  the 
rigid/perfectly-plastic  model  to  predict  Pq.  If  more  detailed 
information  is  required  concerning  deformations  under  a given  load 
P,  use  the  elastic/perfectly-plastic  model  if  P < Pq,  and  use 
the  rigid/strain-hardening  model  if  P > P^. 

In  the  remainder  of  this  paper  we  shall  mostly  be  concerned 
with  the  perfectly-plastic  models,  and  with  their  application 
to  structures  of  various  complexity.  We  shall  begin  with  beams 
and  frames  under  bending  in  Secs.  2-4,  add  tension  in  Sec.  5,  and 
then  consider  torsion  and  torsion  with  tension  in  Secs.  6 and  7. 

We  shall  discuss  both  the  rigid/perfectly-plastic  mdoel  for 
determining  Pq,  and  the  elastic/perfectly-plastic  model  for 
discussing  the  response  at  loads  less  than  Pq.  Sections  8 and  9 
will  then  present  more  general  theories  for  perfectly-plastic 
and  strain-hardening  materials,  respectively,  and  Sec.  10  will 
analyze  some  circular  plates  for  a variety  of  models.  Section  11 
will  present  some  of  the  problems  which  exist  in  applying  finite 
element  methods  to  plasticity  problems,  and  the  paper  will  close 
with  a brief  concluding  section. 

2.  Beams  and  Frames  . We  consider  first  the  pure  bending  of 
an  elastic/perfectly-plastic  beam.  We  make  the  fundamental 
assumption  of  Euler-Bernouilli  beam  theory  that  plane  sections 
remain  plane  and  normal,  whence  the  strain  across  the  section 
will  remain  linear.  For  sufficiently  small  strains  the  stress 
will  also  be  linear.  Fig.*  5a,  but  for  a certain  moment 

* 

Figure  5 and  several  other  figures  are  taken  from  [3]  with 
permission  of  the  publisher.  Such  figures  will  be  identified 
by  an  asterisk  following  their  first  reference. 
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the  stress  will  reach  in  magnitude  at  the  outer  fibers. 

Fig.  5b.  Greater  strains  in  the  outer  fibers  can  still  be 
accomodated  by  |al  = thus  leading  to  the  stress  distribution 
in  Fig.  5c.  In  the  limit,  the  strain  will  become  infinite 
under  the  finite  moment  obtained  from  the  stress  distributon 
in  Fig.  5d.  Figure*  6 shows  the  resulting  moment-curvature 
relation.  Curve  OABD  is  the  curve  for  the  elastic/perf ectly- 
plastic  section  being  considered,  for  the  particular  case  of  a 
rectangular  cross  section.  Curve  OABC  shows  the  diagram  for  the 
mild  steel  of  Fig.  la.  However,  in  the  same  spirit  in  which  the 
nonlinear  curve  in  Fig.  lb  was  replaced  by  the  piecewise  linear 
one,  it  is  usual  to  solve  beam  and  frame  problems  according  to 
OED  or  OFD  depending  on  whether  or  not  elastic  strains  are  to 
be  included. 

The  effect  of  the  "knee"  AB  will  depend  upon  the  shape  of 
the  section.  It  will  be  greater  than  in  Fig.  6 for  a circular 
or  cruciform  section,  less  for  an  I-beam  or  sandwich  section. 
Indeed,  in  the  limit  as  the  bending  resistance  of  web  or  core 
becomes  negligible  even  while  the  flange  or  sheet  thickness 
goes  to  zero,  curve  OED  would  become  exact. 

Figure*  7 shows  a statically  indeterminate  propped  cantilever 
beam  under  increasing  load.  For  F £ the  beam  is  elastic  and 
at  F = Fj^  the  moment  at  A is  “Mq.  As  F is  increased  above  Fj^, 
the  moment  at  A must  remain  at  -Mq , so  that  the  increment  of 
load  above  Fj^  produces  the  same  response  as  if  a hinge  were 
inserted  at  A.  At  F = Fg , M = Mq  at  B.  Since  the  insertion  of 
a second  yield-hinge  at  B would  render  the  beam  a mechanism,  Fg 
is  the  yield-point  load.  Notice  how  closely  this  example  parallels 
the  truss  in  Sec.  1. 
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Now,  if  we  are  interested  only  in  predicting  Fq,  we  may 
proceed  directly  to  the  moment  distribution  of  Fig.  7e.  Indeed, 
since  the  moment  varies  linearly  along  AB  and  BC,  and  since  2 
hinges  are  necessary  to  produce  a mechanism,  it  is  obvious  that 
we  must  have  hinges  at  both  A and  B at  the  yield-point  load. 
Statics,  or  even  more  simply  the  Principle  of  Virtual  Work, then 
leads  immediately  to 

Fq  = 3Mq/L  (2.1) 

A less  trivial  excimple  is  furnished  by  the  portal  frame  in 
Fig.*  8a.  This  frame  can  be  turned  into  various  mechanisms 
by  the  insertion  of  three  yield  hinges  in  four  possible  locations 
as  illustrated  by  Figs.  8b,c,d.  The  corresponding  yield-point 
loads  are  all  different; 

Fj^  = 2Mq/L  F^  = Mq/2L  F^  = 5Mq/8L  (2.2) 

Figures  8e,f,g  show  the  corresponding  moment  distributions. 
Evidently  the  moments  corresponding  to  Fj^  and  F^  exceed  Mq, 
so  that  F^  must  be  the  correct  yield-point  load. 

We  notice  that  F is  the  smallest  of  the  three  mechanism 

c 

loads  so  that  any  of  (2.2)  provide  an  upper  bound  on  the  yield- 
point  load.  This  fact  is  not  just  a coincidence.  There  are 
two  basic  theorems  of  limit  analysis,  proved  independently  by 
Gzodev  [4],  Hill  [5,6],  and  Drucker,  Greenberg,  and  Prager  [7,8,9] 
which  enable  us  to  find  both  upper  and  lower  bounds  on  the  yield- 
point  load. 

These  theorems  are  easily  stated  in  a quite  general  form. 

We  are  given  a structure  with  a load  distribution  fully  defined  by 
a single  parameter  P.  We  define  a statically  admissible  field  as 
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a set  of  moments  which  are  in  internal  and  external  equilibrium 
with  a load  P and  which  nowhere  violate  the  yield  condition.  A 
kinematically  admissible  field  is  defined  by  a mechanism.  At 
every  yield  hinge  a moment  M is  assigned  equal  to  + Mq  in  the 
same  sense  as  the  rotation  6.  Finally  a load  is  defined  so 
that  the  external  work  done  by  P"*"  on  the  mechanism  motion  is 
equal  to  the  internal  work  ZM6 . 

At  the  actual  yield-point  load  P^  all  conditions  of  both 
static  and  kinematic  admissibility  will  be  satisfied.  The 
theorems,  which  are  easily  proved  [3],  state  that  Pq  is  the  largest 
P and  the  smallest  P"*": 

P"  < Pq  1 (2.3) 

xne  importance  of  these  theorems  is  increased  by  the  fact  that 
Pq  in  (2.3)  is  the  same  yield-point  load  for  any  of  the  elastic/ 
perfectly-plastic , rigid/perf ectly-plastic , or  rigid/strain- 
hardening models. 

Examples  of  the  application  of  the  limit  analysis  theorems 
to  rectangular  and  non-rectangular  frames,  Vierendeel  trusses, 
and  grids  may  be  found  in  [3,  Chap.  3]. 

3.  Unloading  and  reloading.  Thus  far  we  have  considered 
only  monotonically  increasing  loads  and  stresses.  If  the  stress 
is  decreased  in  an  elastic  material,  the  stress-strain  curve  will 
simply  be  retraced.  However,  for  an  elastic/perf ectly-plastic 
material  loaded  into  the  plastic  range  to  point  B in  Fig.  9, 
unloading  would  not  retrace  BAO,  but  would  be  along  BC  parallel 
to  the  initial  elastic  line  OA.  Reloading  from  C would  retrace 
CD  and  continue  along  BD  as  if  the  unloading  excursion  had  not 
taken  place. 
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If,  instead  of  reloading  at  C,  unloading  were  to  continue 
until  the  stress  were  zero,  a permanent  strain  OE  would  remain 
in  the  specimen  under  zero  stress.  Further  unloading  would 
continue  along  the  same  line  BCE  until  the  compressive  yield 
stress,  -Oq  was  reached  at  F.  The  strain  could  further  decrease 
to  G with  no  change  in  stress.  Reloading  at  this  point  would  be 
along  GH  parallel  to  BF  and  the  original  elastic  line  OA. 

It  is  clear  from  the  above  description  that  neither  stress  nor 
strain  is  a unique  function  of  the  other,  so  that  the  constitutive 
relations  must  be  expressed  in  rate  form,  with  a need  for  later 
generality  we  may  write  the  law  as  follows: 

a ~ £ = £ +e^  (3.1a,b,c) 

= E a = 2Xa  (3.1d,e) 

IF  a - Cq  <0  OR  2oa < 0 THEN  A = 0 (3. If) 

Observe  that  when  a remains  at  e is  not  determined  but  must 

be  non-negative,  whereas  if  a remains  at  e ^ 0.  For  any 

given  history  of  strain  it  will  always  be  clear  which  branch 
of  (3.1)  is  relevant,  so  that  the  equations  may  be  trivially 
integrated. 

The  rigid/perf ectly-plastic  model  is  a special  case  of  the 
above  in  which  the  elastic  portions  are  vertical  in  Fig.  9.  The 
equations  are  given  by  (3.1)  with  E = “ in  ( 3 . Id) , whence 
= 0 in  (3.1c). 

Unloading  and  reloading  for  strain-hardening  materials  is 
more  complicated  and  we  delay  the  subject  until  Sec.  9. 

Some  of  the  simplified  results  and  techniques  used  in  plastic 
structural  analysis  are  valid  only  in  the  absence  of  plastic 
unloading.  For  this  reason,  it  is  important  to  note  that  local 


unloading  may  take  place  even  under  a monotonically  increasing 
external  load.  A simple  example  of  this  phenomenon,  first 
discussed  by  Drucker  [10] , is  shown  in  Fig.  10.  Three  elastic/ 
perf ectly-plastic  bars,  with  relative  yield  strengths  indicated 
by  the  circled  numbers,  are  attached  to  a fixed  rigid  bar  AB  and 
an  otherwise  free  rigid  bar  CD.  As  P is  increased,  all  bars  will 
first  be  elastic  in  tension  and  bar  1 will  yield  first  at  P = 

Bars  2 and  3 will  have  increasing  tensile  stresses  as  P is 
increased  above  P^^  to  P2.  For  P > P2  the  load  increment  must 
be  carried  by  bars  1 and  2 so  that  a compressive  increment  is 
added  to  bar  1,  i.e.,  it  unloads.  Eventually,  at  P = Pq,  bar 
1 yields  in  compression  and  a mechanism  is  formed.  Thus,  even 
though  P is  monotonically  increased  to  Pq,  the  stress  in  bar  1 
is  first  increased  to  Oq,  remains  there  while  its  strain  increases 
and  then  decreases  to  -Oq  when  failure  occurs.  A simple  portal 
frame  which  exhibits  similar  behavior  is  discussed  in  [3,  Chap.  4] 
and  a more  complicated  truss  with  several  such  stress  reversals 
is  analyzed  in  [11]  . 

4.  Other  limitations . The  yield-point  load  will  always  be 
a limitation  on  the  maximum  load  which  a structure  can  support. 
However,  there  are  circumstances  in  which  other  considerations 
may  be  more  restrictive.  For  example,  slender  members  subject 
to  compressive  load  may  buckle  elastically  or  even  plastically 
before  the  yield-point  load  is  reached. 

Another  possibility  is  excessive  deformations.  Consider, 
for  example  [12]  the  beam  on  four  supports  in  Fig.*  11.  The 
yield-point  load  is  4Mq/L,  independently  of  the  distance  CL 
between  the  outer  supports.  For  very  small  C the  configuration 
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is  similar  to  a clamped  beam  of  length  2L  which  has  the  same 
yield-point  load.  However,  for  large  C the  effective  strength 
added  by  the  outer  supports  would  be  expected  to  decrease.  Indeed, 
in  the  limit  as  ;;  becomes  infinite  the  central  span  should  act 
like  a simply  supported  beam  whose  yield-point  load  is  only 

2Mq/L. 

The  solution  of  this  paradox  lies  in  considering  the 
displacements.  Figure*  12  shows  the  displacements  just  prior 
to  the  yield-point  load  for  various  For  Q = 10  the  displacements 
would  almost  certainly  render  the  beam  unservicable  before  the  load 
has  reached  its  yield-point  load.  For  example,  if  the  maximum 
allowable  displacement  is  5 times  the  maximum  elastic  dis- 
placement in  the  built-in  beam,  the  yield-point  load  is  the 
limiting  criterion  only  if  c,  £ 2.  For  C = 5 the  maximum  dis- 
placement is  reached  when  F = 3Mq/L,  for  i;  = 10  it  occurs  when 
F = 2.5Mq/L.  In  the  limit,  as  C tends  to  infinity,  the  maximum 
allowable  displacement  occurs  with  the  formation  of  the  first 
hinge  at  F = 2Mq/L  in  agreement  with  the  simply  supported  beam. 
Details  of  the  solution  may  be  found  in  [3,  Chap.  4]  together 
with  a simplified  method  for  estimating  displacements. 

A less  obvious  type  of  limitation  can  occur  when  there  is  more 
than  one  load  parameter.  Consider,  for  example,  the  frame  in 
Fig.*  13,  under  independent  loads  F^  and  F2.  The  methods  of 
limit  analysis  are  readily  adapted  to  this  two-parameter  situation, 
and  the  results  are  shown  in  Fig.*  14.  Load  points  f^  = LF^/M^ 
lying  on  the  hexagon  MNPQRS  correspond  to  yield-point  with  each 
side  representing  a different  mechanism,  points  inside  the 
hexagon  do  not  yield,  and  points  outside  are  unattainable. 
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Now  suppose  that  £2  is  increased  to  point  A corresponding 
to  a fixed  dead  weight  load,  and  is  then  alternated  between 
A and  B corresponding  to  a variable  wind  load.  The  plastic 
part  of  the  resultant  displacements,  computed  according  to  an 
approximate  but  complete  elastic-plastic  analysis,  is  shown  in 
Fig.*  15.  Here  6^^  refers  to  the  dimensionless  sidesway, 

62P  to  the  dimensionless  displacement  under  the  load,  and  the 
actual  dimensioned  quantities  are  given  by 

= (MqL^/6EI)5^  (4.1) 

For  the  rigid/plastic  model,  E is  infinite  so  that  = 0 
for  any  number  of  cycles.  However,  for  any  finite  E,  however 
large,  A^  will  become  arbitrarily  large  after  a sufficient 
number  of  cycles,  so  that  the  structure  can  support  only  a 
limited  number  of  cycles  even  though  the  load  is  always 
strictly  within  the  yield-point  hexagon  of  Fig.  14. 

Even  one  load  can  cause  undesirable  behavior  if  it  is 
alternated.  Suppose,  for  example,  that  there  is  no  f^,  but 
that  f2  alternates  between  points  A and  C.  The  resulting  plastic 
displacement  is  shown  by  the  dashed  curve  in  Fig.  15,  and  is 
not  in  itself  alarming  since  it  remains  finite.  However, 
each  increase  in  6^^  corresponds  to  a positive  plastic  rotation 
of  the  yield  hinge  under  the  load,  whereas  each  decrease  is  a 
negative  plastic  rotation.  Construction  of  a moment-rotation 
diagram  would  show  that  a large  hysterisis  loop  is  present, 
representing  substantial  energy  input  to  the  hinge.  Although 
details  of  the  phenomenon  cannot  be  discjssed  within  any  of  the 
macroscopic  theories  considered  here,  it  is  evident  that  a 
relatively  few  cycles  may  cause  internal  damage  which  would 
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render  the  frame  unserviceable . 
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The  shakedown  load  of  a structure  is  defined  as  the  largest 
load  for  which  none  of  the  undesirable  features  described  above 
occurs.  Techniques  are  available  for  finding  or  estimating  the 
incremental  collapse  load  and  the  cyclic  collapse  load  [3,  Chap. 
The  smallest  of  these  two  and  the  yield-point  load  is  the 
shakedown  load. 


5] 


Beams  under  tension  and  bending.  If  a beam  is  subject 


to  both  axial  force  and  bending  moments,  the  strain  rate  at 
any  distance  z from  i-he  center  axis  is 

£=e+z<  (5.1) 

Elastic-plastic  stages  will  be  similar  to  Fig.  5,  except 
that  the  neutral  and  center  axes  will  generally  be  different. 

For  a rigid/perfectly-plastic  material  we  are  interested  only 
in  a fully  plastic  cross  section.  Here,  there  are  four 
possibilities  as  shown  in  Fig.  16.  The  corresponding  resultants 


for  a rectangular  section  must  satisfy 


N = -Nq^ 

M = Mq(1-^^  ) 

• 

< 

> 

0 

1 

A 

A 

1 

(5.2a) 

N = Ng? 

M = -Mq(1-^^) 

• 

IC 

< 

0 

-1  < ^ < 

1 

(5  .2b) 

M = -Ng 

M = 0 

• 

e 

< 

0 

Ui  1 1 

(5.2c) 

N = Nq 

M = 0 

• 

e 

> 

0 

Ul  i 1 

(5. 2d) 

where  ?h  is  the  distance  to  the  neutral  axis.  In  addition, 


e and  < must  be  related  by 

e+c<=0  (5.3) 

since  the  neutral  axis  strain  rate  is  zero  by  definition. 

Equations  (5.2)  define  an  interaction  curve  in  stress- 
resultant  space,  as  shown  in  Fig.  17.  Points  on  the  curve 
correspond  to  a fully-plastic  section,  points  within  correspond 


to  a section  at  least  partly  elastic,  and  points  outside  the 
curve  are  unobtainable. 

Let  Q = (N,M)  denote  a generalized  stress  vector,  and 
q = (e,K)  a general  strain-rate  vector.  The  internal  dissipation 
rate  of  the  beam  can  then  be  written 


D 

int 


L , . L ^ 

f aedv  = / (Ne  + MK)dx  = f Q*q  dx 
V 0 0 ~ 


(5.4) 


so  that  Prager's  criterion  [13]  for  genralized  variables 
satisfied. 

Now,  except  for  A and  D,  each  point  Q on  the  yield  curve 
corresponds  to  a unique  <;  and  hence,  according  to  (5.3)  to  a 
unique  direction  for  q.  Further,  it  is  easily  shown  that  q 
is  perpendicular  to  the  yield  curve.  At  points  A and  D 
neither  ^ nor  the  normal  to  the  curve  are  unique.  However  al^ 
admissible  values  of  ? correspond  to  directions  between  the 
limiting  normals  so  that  in  a sense  we  may  state  that  the 
generalized  strain-rate  vector  is  always  normal  to  the  yield 
curve.  The  length  of  q is,  of  course,  undetermined. 

To  express  the  above  geometric  results  algebraically,  we 
eliminate  ^ between  (6.2a)  or  (6.2b)  to  obtain  the  equations 
for  arcs  DBA  and  DEA,  respectively: 


fj^(N,M)  = (N/Nq)^  + M/Mq  -1  = 0 

f2(N,M)  = (N/Nq)^  - M/Mg  -1=0 

If  q is  normal  to  (5.5a)  then 

e = X9f/3N  = 2AN/Nq^ 

K = Xaf/9M  = VMg 

where  X is  non-negative.  In  vector  form. 


(5.5a) 

(5.5b) 

(5.6) 

the  strain  rates  for 


DBA  and  DEA  are,  respectively 
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where  a and  u are  each  non-negative.  Finally,  at  points  D or 
A the  requirement  that  q lie  between  the  limiting  normals  can 
be  expressed  by 

q = XVf^  + )jyf2  A > 0 y > 0 (5.8) 

The  above  theory  may  be  used  to  assess  the  influence  of 
axial  forces  on  the  yield-point  loads  of  frames,  or  to 
analyze  problems  of  curved  beams  or  arches.  Some  applications 
may  be  found  in  [3,  Chap.  7]. 

6.  Torsion  of  a circular  cylindrical  bar.  A circular 
cylinder  of  length  L is  fixed  at  one  end  and  subjected  to  a 
pure  torque.  We  make  the  usual  assumption  that  each  cross 
section  undergoes  only  a rotation  proportional  to  its  distance 
from  the  fixed  end.  Then  the  only  displacement  is  the 
tangential  one 

V = azr  (6.1) 

where  a is  the  angle  of  twist  per  unit  length.  The  only  non-zero 
strain  is  the  shearing  strain 

Y = Yqz  = (6.2) 

For  a sufficiently  small,  the  entire  bar  will  be  elastic 

with 

T = Gar  (6.3) 

The  resultant  torque  in  this  case  is 

a . 

T = 27t  / rtr  dr  = (TT/2)Gaa^ 

0 


(6.4) 
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As  a is  slowly  increased,  t will  reach  its  yield  value  k 
in  shear  at  the  boundary  of  the  bar  for  a value  defined  by 


G a = k 

so  that  the  maximum  elastic  torque  is 

= (TT/2)ka^ 


(6.5) 


(6.6) 


For  larger  values  of  a , an  outer  annulus  of  the  bar  will 
become  plastic  while  the  inner  core  remains  elastic.  Thus, 


for  a perf ectly-plastic  material 


T=Gar  0<r<p 


(6.7) 


T = k 


p < r < a 


Since  x must  be  continuous  at  the  elastic-plastic  boundary  p , 


P = k/Ga 


The  torque  in  this  case  is 


T = I .ica^U  - i (^)^] 


Evidently  the  limiting  fully  plastic  torque 

Tq  = (2/3)  TTka^ 


(6.8) 


(6.9) 


(6.10) 


is  attained  only  as  a tends  to  infinity.  However,  we  note 
that  T is  within  less  than  1 percent  of  Tq  when  a is  only  three 
times  its  maximum  elastic  value 


7.  Combined  tension  and  torsion.  If  a beam  of  circular 
cross  section  is  subjected  to  both  tension  and  torsion,  then 
it  will  certainly  be  limited  by  the  tensile  and  shear  yield 
stresses  -Jq  and  k.  These  values  are  both  experimentally 
determined  constants  so  that  one  can  be  expressed  in  terms  of 
the  other,  and  there  is  considerable  evidence  that  =/3c . 
Thus  the  stress  point  Q = (a,x)  must  certainly  lia  within  the 
rectangle  ABCD  in  Fig.  18.  However,  if  a and  i are  both  non- 
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zero,  experience  suggests  that  their  interaction  will  pose  a 
more  severe  limitation  and  that  Q must  lie  within  some  curve 
within  the  rectangle.  In  the  next  section  we  will  justify 
taking  this  curve  to  be  the  ellipse 

f(Q)  = [a^/2  + -1=0  (7.1) 

Let  us  visualize  an  experiment  in  which  the  unit  extension 
e and  unit  rotation  a are  controlled.  In  terms  of  these 
parameters,  the  strain  vector  is 

q s (e,Y)  = (e,ar)  (7.2) 

For  e and  a sufficiently  small,  the  bar  will  be  everywhere 
elastic,  hence  the  stress  vector  is 

Q E (a,T)  = G(3e,ar)  (7.3) 

where  we  have,  for  simplicity,  taken  v = 1/2  so  that  E = 3G. 
Equation  (7.3)  will  be  valid  so  long  as  Q lies  everywhere 
within  the  ellipse  (7.1),  i.e.,  so  long  as 

(/3Ge/k)^  + (Gaa/k)^  - 1 < 0 (7.4) 

If  e and  a are  such  that  part  of  the  bar  is  plastic,  the 
constitutive  equations  must  be  expressed  in  terms  of  rates. 

To  this  end,  we  first  rewrite  Eqs.  (3.1)  in  vector  form  and 


apply  them  to  the  present  problem: 

2 

f(Q)  s (1/k  ) (a^/3  + T^)  - 1 < 0 (7.5a) 

X > 0 q = q®  + (7.5b,c) 

q®  = G'^(o/3,t)  (7.5d) 

qP  = XVf  = 2(X/k^) (a/3,T)  (7.5e) 

IF  f < 0 OR  f < 0 THEN  X = 0 (7.5f) 


In  generalizing  (3.1e)  we  have  been  guided  by  the  results  of 


Sec.  5. 


Let  us  consider  specifically  the  case  where  X > 0.  It 
then  follows  from  (7.5f)  that  f = f = 0,  whence  the  local 
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dissipation  rate  is 

Q • q = + TT)  + + T^)  = 2X  (7.6) 

We  can  now  rewrite  Eqs.  (7.5)  in  the  more  useable  form 

f (Q)  ^ 0 Q*q  ^ 0 (7.7a,b) 

IF  f < 0 OR  f < 0 THEN 


ELSE 


0 = 3Ge  T = Gra  (7.7c) 

0 = (G/k^)[3k^  - 0^)e  - 0ira] 

(7.7d) 

T = (G/k^) [-0Te  + (k^  - T^)ra] 


Let  us  consider  a specific  loading  sequence  in  which  e 
is  slowly  increased  from  0 to  k//3G  and  is  then  held  there 
while  a is  increased  to  k/Ga.  The  first  stage  is  fully 
elastic  and  is  governed  by  (7.7c)  or,  since  we  start  from 
zero,  by  (7.3).  At  the  end  of  this  stage 

e = k//3G  0 = kv'l’  a = T = 0 (7.8) 


At  this  point  f = 0 for  all  r,  hence  (7.7d)  are  to  be 
integrated  with  (7.8)  as  initial  conditions.  This  process 
leads  to 

0 = k/T  sech  (Gar/k)  x =k  tanh  (Gar/k)  (7.9) 

In  particular,  the  final  stress  state  is 

0 = k/T  sech  (r/a)  t =k  tanh  (r/a)  (7.10) 

Consider  now  an  alternative  loading  sequence  which  reaches 
the  Scune  final  state  by  first  increasing  a to  its  final  value, 
and  then  increasing  e.  As  shown  in  [14,  Sec.  14],  the  final 
stress  state  in  this  case  is 


a ^ /I  tanh  ^ x =k  sech  ^ 

C = 1 - u + tanh“^  u u = /I  - r^/a* 


(7.11) 


8.  Multi-axial  stress  states.  The  most  general  state  of 
stress  is  characterized  by  the  six  components  of  the  symmetric 
tensor  a.  Therefore,  general  behavior  of  a perfectly-plastic 
material  will  be  characterized  by  a yield  condition 

f(a)  < 0 (8.1) 

If  the  material  is  isotropic,  the  form  of  the  yield  condition 
cannot  depend  upon  the  orientation  of  the  coordinate  axes,  so 
that  f can  have  only  three  independent  arguments.  Thus  (8.1) 
can  be  replaced  by  either  of 

f(c,  J2,  J^)  ± 0 f(a^^,  02t  a^)  £ 0 (8.2a,b) 
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where  a is  the  mean  normal  stress  and  and  are  the  second 

and  third  invariants  of  the  stress  deviation  tensor,  and  a. 

1 

are  the  three  principal  stresses.  Experience  indicates  that 
the  mean  normal  stress  does  not  affect  the  yielding  of  most 
structural  materials  whence  (8.2  ) can  be  further  simplified  to 
f(J2»  J3)  1 0 f(a^-C2,  ° (8.3a,b) 

In  1868,  Tresca  [15]  suggested  that  yielding  was  governed 
by  the  maximum  shearing  stress,  hence 

f = max  - 02!,  - a^l,  | ] - 2k(8.4) 

This  condition  is  generally  simple  to  apply  when  the  principal 
stress  directions  are  fixed  and  known  but  is  difficult  in  more 
general  cases.  Therefore,  in  1913,  Mises  [16]  proposed  approximating 
(8.4)  by  a special  case  of  (8.3a): 


f * 6JJ  - 2a„ 

-<“x  - '“x  - ^ <“y 

+ 6(t  ^ + t ^ + t - 2a  ^ 

xy  xz  yz  ) 0 


(8.5) 


If  Oq  = /3k,  then  the  two  yield  conditions  will  agree  in  pure 
shear  and  will  have  their  maxiirium  discrepancy  of  2//3  = 1.155 
in  simple  tension. 

Later  experimental  evidence  has  shown,  that  Mises ' 
yield  condition  is  a closer  approximation  than  Tresca's  to 
most  real  structural  materials,  but  in  view  of  the  small 
difference  between  them,  it  is  customary  to  assume  whichever 
one  is  simplest  in  a given  problem  situation. 

We  have  already  derived  the  constitutive  equations  for  a 
special  case  in  Sec.  5,  and  have  assiimed  them  in  another 
special  case  in  Sec.  7.  In  general  they  are 


f(a)  <0  A>0  e = e +£^ 


(8 . 6a ,b,c) 


e®  = E a eP  = XVf 


(8.6d,e) 


IF  f < 0 OR  f < 0 THEN  X = 0 


(8.6f) 


For  the  particular  case  of  the  Mises  yield  condition  (8.5), 
the  general  results  may  be  simply  expressed  in  terms  of  the 
stress  and  strain  deviators 


s = a - Ic 


e = e - le 


(8.7) 


where  a and  e are  the  mean  normal  stress  and  strain: 


f = s:s  -2k'"<0  X>0 


(8.8a,b) 


e - s/2G  + 


a = 3K£  ' 


(8.8c,d) 
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together  with  (8.6f).  For  X > 0,  the  requirement  f = f = 0 shows 
that 


e;s  = s:s/2G  + Xs:s  = 2k  X 


(8.9) 


whence  (8.8)  can  be  rewritten 


IF  f < 0 OR  f < 0 THEN  s = 2Ge 


ELSE  s = 2G(e  - s(e:s)/2k' 


(8.10) 


Note  that  although  for  a perfectly  plastic  material  the  strain 
rate  is  not  uniquely  determined  by  the  stress  rate,  as  shown 
by  (8.10)  a given  strain  rate  does  uniquely  determine  the 
stress  rate. 

Equations  (8.8),  without  the  elastic  terms,  were  first 
used  by  Mises  in  1913  [16] , following  ideas  advanced  by 
St.  Venant  [17]  and  Levy  [18]  in  1870.  The  elastic  terms 
were  incorporated  by  Prandtl  [19]  and  Reus  [20]  in  1924  and 
1930.  The  more  general  (8.6)  were  first  proposed  by  Mises  in 
1928  [21] . Later,  additional  reasons  for  using  them  were 
advanced  by  Hill  [22]  , Drucker  [23] . 

For  structural  problems  which  are  coveniently  handled  with 
generalized  variables,  it  is  merely  necessary  to  replace  the 
stress  matrix  o by  the  generalized  stress  vector  Q,  and  the 
strain  e by  the  generalized  strain  q.  The  appropriate  form  of 
the  yield  condition  f and  the  elastic  matrix  E must,  of  course, 
be  determined  for  the  particular  structure. 
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9.  Strain-hardening  material.  If  a strain-hardening 
material  is  loaded  to  a point  B in  the  plastic  range.  Fig.  20, 
and  then  unloaded,  the  unloading  will  take  place  along  a line 
BC  parallel  to  the  original  elastic  line  OA.  If  the  load  is 
again  reversed  at  C,  CB  will  be  retraced  and  further  loading 
will  take  place  along  BD  just  as  if  the  excursion  BCB  had  not 
occured.  Actually,  BCB  will  constitute  a small  hysteresis 
loop  and  the  reloading  may  start  to  behave  plastically 
shortly  before  B and  rejoin  line  AD  shortly  after  B,  but  these 
effects  are  small  compared  to  the  overall  elastic  and  plastic 
effects  and  hence  are  ignored  in  most  models. 

If  the  bar  is  fully  unloaded  to  point  E,  a residual  strain 
OE  will  remain  and,  with  further  unloading,  plastic  flow  in 
compression  will  begin  at  a point  F and  continue  along  a line  FG 

Thus  far  the  description  of  behavior  is  an  obvious 
generalization  of  Fig.  9 for  a perf ectly-plastic  material, 
but  two  questions  remain  to  be  answered:  where  is  point  F 

and  what  is  the  slope  FG?  If -elastic  slopes  and  initial 
yield-stress  are  originally  symmetric  in  tension  and  compression 
it  seems  reasonable  to  assvime  that  FG  is  parallel  to  AB.  With 
this  assumption,  the  future  behavior  of  a tensile  specimen 
at  a generic  point  C and  with  a given  loading  history  is 
specified  by  four  numbers:  the  elastic  slope  (OA,  FB) , the 
plastic  slope  (AB,  GF) , the  current  yield  stress  in  tension 
(ordinate  of  B) , and  the  current  yield  stress  in  compression 
(ordinate  of  F) . The  two  slopes  are  material  constants 
independent  of  history,  and,  for  the  history  shown  in  Fig.  20, 
the  current  tensile  yield  stress  is  the  value  at  which 
the  most  recent  plastic  flow  was  reversed.  However,  various 
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models  have  been  proposed  to  predict  the  ordinate  of  F. 

Historically,  most  of  the  early  work  in  plasticity  which 
accounted  for  strain-hardening  at  all  assumed  that  it  was 
isotropic , i.e.,  that  the  current  yield  stresses  in  tension 
and  compression  were  always  equal.  Thus,  in  Fig.  20,  compressive 
plastic  flow  would  take  place  along  • 

However,  experimental  evidence  generally  suggests  that 
strain  hardening  in  tension  will  reduce,  rather  than  raise, 
the  magnitude  of  the  compressive  yield  stress.  To  account 
for  this  "Bauschinger  effect,"  Prager  [24,  25]  proposed  a theory 
which  has  become  known  as  kinematic  hardening.  In  terms  of 
the  tension  test,  kinematic  hardening  assiames  that  the  elastic 
range  remains  constant  at  2Cq,  independently  of  a^.  Thus 
unloading  would  take  place  along  in  Fig.  20.  We  note  that 

if  the  hardening  is  so  great  that  > 2a^,  then 
plastic  flow  in  compression  would  take  place  under  an  unloading 
stress  which  was  still  positive. 

These  two  models  are  generally  accepted  as  the  extreme 
possibilities,  and  any  number  of  intermediate  models  can  be 
proposed.  For  example,  F2G2  in  Fig.  20  shows  a compressive 
yield  stress  which  is  independent  of 

For  multi-axial  status  of  stress/  Eqs.  (8.6)  may  still 
be  used  with  some  modifications  in  the  interpretation  of  f(c) 
and  X.  The  yield  function  f(a)  is  no  longer  a fixed  property 
of  the  material  but  will  vary  with  stress  history.  It  will 
remain  constant  when  the  behavior  is  elastic,  but  will  vary 
during  plastic  flow  so  that  a is  always  on  the  yield  surface. 

Thus  one  should  write  f = f(a,H)  where  H stands  symbolically 
for  the  entire  stress-strain  history.  Then,  in  Eq.  (8.6f), 
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f is  not  a total  derivative  of  the  entire  function,  but  is 
defined  specifically  by  f = 7f  • a. 

Since  the  requirement  that  f = 0 during  plastic  flow  is 
a partial  constraint  on  the  history  of  f rather  than  the 
equation  f = 0 , we  have  one  fewer  constitutive  equation  available. 
The  system  is  balanced  by  requiring  a relation  between  the 
plastic  strain-rate  magnitude  X and  the  yield  function 
increase  f.  Thus,  when  X is  not  equal  to  zero,  it  is  given 
by 

X = Bf  (9.1) 

For  a hardening  curve  with  variable  slope,  B could  be  a function 
of  stress;  for  the  linear  hardening  here  considered,  B is  a 
material  constant  proportional  to  the  slope  of  AB  in  Fig.  20. 

There  still  remains  the  problem  of  how  f varies  with 
plastic  loading.  One  reason  for  the  popularity  of  the  isotropic 
and  kinematic  models  is  that  they  both  give  easily  implemented 
answers  to  that  question.  According  to  isotropic  hardening, 
the  yield  surface  remains  fixed  in  origin,  orientation,  and 
shape,  but  will  change  in  size  to  accomodate  plastic  flow.  Thus, 
if  the  initial  yield  function  is  written 

f(j)  = g(a)  - Cq  (9.2) 

the  current  yield  function  after  plastic  flow  will  be 

f(a,  H)  = g(a)  - (9.3) 

where  is  the  largest  value  of  g(a)  attained  during  the 
preceding  history. 

In  geometric  terms , kinematic  hardening  states  that  the 
yield  function  remains  fixed  in  orientation,  shape,  and  size, 
but  will  change  its  origin  to  accomodate  plastic  flow.  Thus 
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Besseling  [27]  and  applied  by  Vaugh  [28]  . Each  material  element, 
infinitesimal  or  finite,  is  replaced  by  a number  of  subelements 
connected  in  parallel-  If  each  sub-element  is  elastic/perfectly- 
plastic  but  with  different  yield  stresses,  the  resulting  total 
model  has  many  similarities  to  kinematic  hardening,  except  that 
the  resulting  B is  piecewise  constant  rather  than  a single 
number.  By  including  some  subelements  with  isotropic  hardening, 
models  may  be  constructed  which  are  intermediate  to  overall 
isotropic  or  kinematic  hardening. 

The  subject  of  proper  strain-hardening  models  is  one  of 
intense  current  interest.  Indeed,  an  entire  symposium  was 
devoted  to  the  topic  at  the  1976  ASME  meeting  and  the  interested 
reader  is  referred  to  the  Proceedings  of  that  meeting  [29]  for 
some  of  the  current  research  being  done  in  the  area. 


10.  Circular  plate.  For  a circular  plate  under  rotationally 
symmetric  loading,  appropriate  generalized  stresses  are  the 
dimensionless  moments 


(m^,  mg)  = 


H 


-H 


Cg ) zdz 


(10.1) 
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where  = OqH  is  the  yield  moment  and  2H  the  plate  thickness. 
Since  the  principal  directions  are  known  and  fixed,  it  is 
convenient  to  use  Tresca ' s yield  condition.  Setting  the  third 
principal  direction  = 0 in  (8.4)  and  using  (10.1)  we  obtain 

f = max  [jm^l,  |m^|,  |m^  - mg|]  - 1 (10.2) 

as  illustrated  in  Fig.*  21. 

For  a constant  load  the  moments  must  satisfy  the  equilibrium 
equation 

(rm^) " - mg  = - 3pr^  (10.3) 

Evidently  if  a rigid-plastic  plate  is  plastic  corresponding  to 
one  of  the  six  sides  of  (10.2)  , the  problem  is  statically 
determinate  and  easily  solved.  For  example,  the  stress 
profile  for  a simply  supported  plate  of  unit  radius  lies 
everywhere  on  side  BC 

nig  = 1 (10.4a) 

Substitution  of  (lCL4a)  in  (10. 3)  and  use  of  the  boundary 

condition  m = m , at  r = 0 leads  to 
r ^ 

mj.  = 1 - pr^  (10.4b) 

and  the  boundary  condition  m^(l)  = 0 furnishes  the  yield-point 
load 

p = 1 (10.4c) 

Before  accepting  (10.4)  as  the  solution,  two  things  must  be 
verified.  First,  one  must  show  that  the  stress  profile  lies  on 
the  finite  side  BC,  which  is  more  restrictive  than  the  infinite 
line  (10.4a).  Since  (10.4b,c)  show  that  0 ^ m^  £ 1 for  all 
0 <_  r £ 1,  this  condition  is  satisfied. 

Second,  we  must  associate  a kinematically  admissible 


i 

i 


velocity  field  with  (10.4).  Generalized  strains  for  this  problem 
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are 

K = - v''  K.  = - v'/r  (10.5) 

* I r a 

On  side  BC  the  strain-rate  vector  must  be  vertical  and  directed 

upwards . Thus 

v"  = 0 w"  = -C  w = C(1  - r)  (10.6) 

where  C is  any  positive  constant. 

At  the  center  of  the  plate  the  slope  jumps  discontinuously 
from  the  symmetry  value  of  0 to  the  negative  value  -C.  Thus, 
it  follows  from  (10.5a)  that  is  positively  infinite.  Since 
r = 0 is  at  the  corner  B in  Fig.  21,  a positive  is  admissible 
Thus  (10.6)  is  a kinematically  admissible  field  and  hence 
(10.4)  is  a complete  moment  solution. 

A clamped  plate  is  less  simple,  since  the  stress  profile 
lies  on  BC  near  the  center  of  the  plate  and  on  CD  near  the  edge. 
The  solution  in  this  case  is  [3,  Chap.  10] 

2 

0 <_  r : mg  = 1 m^  = 1 - pr  w = Cn(l~log  n-r/^) 

2 

n < r < 1:  mg  = log  r + (3p/2) (1-r  ) m^  = m^  - 1 (10.7a) 

w = - C n log  r 

where  p and  n are  defined  by 

3p  - log  p - 5 = 0 pn^  = 1 (10.7b) 

It  is  readily  verified  that  the  stress  profile  lies  on  the 
finite  sides  and  that  the  velocity  field  is  kinematically 
admissible . 

As  a final  example  of  a rigid/perf ectly-plastic  plate  we 
consider  the  annulus  in  Fig.  22  [30]  . The  stress  profile  in 
Fig.  23a  is  statically  admissible,  but  the  associated  velocity 
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field  shown  in  Fig.  23d  is  not  kinematically  admissible, 
since  <g  < 0 for  1 < r < C-  The  velocity  profile  in  Fig.  23e 
would  be  kinematically  admissible  if  associated  with  the 
stress  profile  in  Fig.  23b,  but  it  turns  out  that  this  profile 
is  not  statically  admissible.  The  actual  solution  is  shown 
in  Fig.  23c, f,  and  involves  an  annulus  of  the  plate  remaining 
rigid.  Details  of  this  solution  along  with  other  rigid/perf ectly- 
plastic  plate  problems  may  be  found  in  [3,  Chap.  10]  and  [31, 

Chap.  4] . 

We  consider  next  an  elastic/perf ectly-plastic  simply 
supported  plate,  using  the  idealized  curve  OED  in  Fig.  6. 

This  is  equivalent  to  considering  an  ideal  sandwich  plate. 

The  general  elastic  solution  for v = 1/3  is  [32] 

h^  nw'  = 4B/3r  + 2Cr/3  + (p/3)r^ 

m = B/r^  - C - (5p/4)r^ 

^ (10.8) 

mQ  = B/r^  - C - (3p/4)r^ 
where  h = H/A,  n = E/Oq,  and  w = W/H. 

For  p sufficiently  small  the  plate  is  everywhere  elastic 


with 

h^  nw"  = (p/6) (-5r  + 2r^) 

, P (10.9) 

ntj.  = (5p/4)  (1  - r^)  m^  = (p/4)  (5-3r^) 

This  solution  remains  valid  until  p = 0.8  when  the  plate  center 
reaches  point  B in  Fig.  21.  For  p slightly  larger  than  0.8, 
the  center  of  the  plate  will  become  plastic  on  side  BC,  with 
the  outer  annulus  remaining  elastic.  Thus  (10.4a,b)  hold 
for  0 ^ r £ ^ and  (10.8)  for  C ^ !•  Using  continuity  at 

r = C we  obtain  the  complete  moment  solution 


y 
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0<r<C:  m^=l  + (p/4)  (-c'^/r^  + 25^  - 5r^) 

mg  = 1 + (p/4)  + 2C^  - 3r^)  (10.10) 

2 

5^r^l:  m^=l-pr  mg=l 

This  solution  holds  until  the  entire  profile  lies  on  BC  at 
the  yield-point  load  p = 1. 

To  find  the  slope  for  p > 0.8,  we  observe  that  in  the 
plastic  regime  BC,  must  vanish,  so  that 

nHK  =nHK®=m  - m„ /3  (10.11) 

r r r 9 

But  (10.11)  also  held  when  the  point  in  question  was  still 
elastic,  so  that  it  can  be  trivially  integrated  with  respect 
to  time  to  yield 


h^  n w"  = - (m^  - mc,/3)  (10.12) 

Therefore,  the  correct  slope  solution  for  p > 0.8  is 

0 ^ r £ r;  h^nw"  = 2r/3  + (p/3)  (-s'^/r  - 5^r  + r^) 

P . (10.13) 

5 <_  r £ 1 h nw'  = 2r/3  + (p/3)  (-2C^  + r"^) 


The  displacement,  of  course,  can  be  found  by  integrating  the 
slope  with  the  boundary  condition  w(l)  = 0.  Curve  (d)  in 
Fig.  24  shows  the  resulting  load  vs.  center  deflection  for 
the  elastic/perf ectly-plastic  plate.  Details  of  the  solution 
may  be  found  in  [32]  . 

This  problem  can  also  be  solved  for  a elastic/strain- 
hardening material.  The  solution  for  p < 0.8  is,  of  course, 
the  fully  elastic  one,  independent  of  any  potential  strain- 
hardening. However,  for  p > 0.8,  two  plastic  regions  begin  to 
form  simultaneously  at  the  plate  center  with  0 £ r £ ^ in  the 
corner  regime  B,  C < r < i;  on  the  side  BC,  and  s r £ 1 


still  elastic.  For  isotropic  hardening  the  equations  can 
be  solved  for  each  regime  and  the  constants  of  integration 
determined  by  continuity  and  boundary  conditions.  For  a value 
of  p slightly  greater  than  1 (the  exact  amount  will  depend 
upon  the  strain-hardening  constant)  , c,  will  equal  one.  For 
larger  p,  the  entire  plate  will  be  plastic  with  0 £ r ^ ? 
in  corner  B and  C £ r ^ 1 on  side  BC.  Algebraic  details  may 
be  found  in  [32  ] . 

Figure  24  shows  the  load-deflection  curves  for  both  the 
elastic/strain-hardening  and  rigid/strain-hardening  models. 

Aside  from  the  fact  that  the  "knee"  in  the  curves  which 
represents  a partially  plastic  plate  is  curved  rather  than 
linear.  Fig.  24  is  qualitatively  very  similar  to  Fig.  4. 
Therefore,  the  same  conclusions  can  be  drawn  from  it,  and  we 
can  gain  confidence  in  our  basic  approach  to  plastic  structures; 
namely  that  the  rigid/perfectly-plastic  model  can  be  used  to 
find  the  yield-point  load,  that  strain-hardening  can  be  neglected 
below  the  yield-point  load,  and  that  elastic  strains  can  be 
ignored  above  it. 

11.  Finite  element  model.  A finite  element  model  has 
become  the  accepted  way  of  solving  complicated  problems  in 
elasticity,  and  usually  the  only  question  is  which  finite 
element  model  to  use.  Essentially,  a trade  off  must  be  made 
between  taking  many  simple  elements  with  simple  internal  fields 
as  opposed  to  fewer  complex  elements  with  complex  fields. 


Almost  any  finite  element  model  for  elasticity  can  be 
easily  adapted  to  include  plasticity  with  or  without  strain- 
hardening. However,  there  are  certain  features  of  plasticity 
which  are  basically  different  from  elasticity  and  which  may 
well  affect  the  desirability  of  any  given  model. 

Since  plasticity  must  be  formulated  as  a rate  problem, 
the  basic  approach  is  as  follows.  At  a generic  time  t^ , it 
is  assumed  that  stress,  strain,  displacement,  and  any  other 
quantities  of  interest  are  all  known,  and  a finite  element 
model  is  used  to  find  the  rates  at  time  t^.  With  all  rates 
known,  time  is  given  a finite  increment  At  and  quantities  at 
time  tg  + At  are  approximated  by 

(j>(tg  + At)  = <j>(tg)+  At  ^ (tg)  (11.1) 

The  process  is  then  repeated  with  tg  replaced  by  tg  + At. 

There  are,  then,  two  distinct  repeated  stages  to  the 
solution:  solving  the  rate  problem  and  applying  (11.1).  We 

consider  first  three  aspects  of  the  rate  problem.  For  definiteness 
we  consider  an  elastic/perfectly-plastic  material  whose 
constitutive  behavior  is  governed  by  Eqs.  (8.6).  The  primary 
problem  is  involved  with  the  branch  points  in  (8.6f)  and  the 
inequalities  (8.6a,b). 

One  of  the  simplest  and  most  useful  finite  element  models 
is  the  constant  strain  one  which  can  be  associated  with 
arbitrary  triangular  or  tetrahedral  elements  in  two  or  three 
dimensions,  respectively.  With  this  element,  f and  A will 
each  have  a constant  value  over  the  entire  element,  so  that 
the  decision  concerning  branch  points  must  be  made  only  on 
an  element-by-element  basis. 


Now,  at  time  tp,  f is  known  so  that  the  first  branch 
point  in  (8.6f)  is  no  problem,  and  when  f < 0 the  element 
is  elastic.  However,  f will  be  obtained  as  part  of  the  rate 
solution  so  that  it  is  not  clear  which  branch  should  be  used 
when  f = 0 . 

A method  which  appears  to  work  well  in  practice,  at  least 
for  one-paraimeter  loads,  is  to  make  any  reasonable  guess 
as  to  whether  each  element  with  f = 0 is  elastic  or  plastic 
and  solve  the  rate  problem.  Each  element  is  then  tested.  If 
it  was  assumed  elastic  with  X = 0,  then  f must  be  non-positive 

or  else  (8.6a)  would  be  violated;  if  it  was  assumed  plastic 

• • 

with  f = 0,  then  X must  be  non-negative  or  else  (8.6b)  would  be 
violated.  Each  element  which  fails  its  test  has  its  assumption 
reversed  and  a new  rate  solution  is  found  and  tested.  This 
process  is  continued  until  all  elements  simultaneously  pass 
their  test.  A simple  example  has  been  constructed  [11]  in  which 
all  tests  are  not  passed  until  the  third  iteration.  Also,  to 
the  author ' s knowledge  there  is  no  proof  that  the  iteration 
process  will  always  converge,  although  there  are  no  known 
counter  examples.  Further  research  is  needed  in  this  area,  and 
meanwhile  it  is  important  that  any  finite  element  program 
continues  to  test  all  assumptions  until  a correct  solution  is 
obtained. 

An  alternative  approach  to  the  rate  problem  is  to  use  the 
plastic  minimum  principles  [33]  and  a mamerical  minimizing 
technique.  The  problem  formulates  as  one  in  quadratic  programming 
[34,  35].  The  question  of  branch  points  is  transformed  to  one 
of  asking  which  constraints  are  active  which  is  automatically 
handled  by  any  quadratic  programming  program. 


The  second  aspect  of  the  rate  problem  lies  in  the  choice 
of  element.  If  any  element  has  variable  stress,  then  at  time 
tg  part  of  the  element  has  f < 0 and  part  of  it  f = 0.  For 
the  f = 0 part,  the  branch-point  questions  raised  above  will 
exist  with  the  further  complication  that  the  correct  rate 
solution  may  be  f = 0 in  part  of  that  domain  and  A = 0 in  the 
remainder.  Even  if  this  problem  is  solved,  numerical 
integrations  across  the  element  will  be  necessary,  thus 
substantially  increasing  the  computer  time  spent  in  finding 
each  rate  solution.  Therefore,  questions  concerning  many 
simple  versus  few  complex  elements  should  be  considered 
directly  for  elastic-plastic  problems,  with  a strong  presumption 
that  simpler  elements  will  be  optimum  than  for  purely  elastic 
structures . 

The  final  aspect  of  the  rate  problem  is  that  of  continuity. 
Elasticity  problems  are  extremely  continuous  and  efficient 
finite  element  methods  will  take  advantage  of  that  fact. 

However,  problems  for  elastic/perfectly-plastic  structures  may 
have  discontinuities  in  stress  derivatives  (as  across  the 
characterisitics  in  a plane  strain  problem)  and  even  of 
stresses  themselves  as  in  Figs.  5d  or  16.  Stress  discontinuities 
pose  no  particular  problem,  the  finite  element  model  merely 
forcing  them  to  occur  only  at  element  boundaries,  but  velocity 
discontinuities  may  also  occur.  For  example  in  Secs.  2 and  10 
we  saw  that  slope  discontinuities  could  occur  in  beams  and 
plates.  It  has  been  shown  [36]  that  finite  element  models 
which  allow  for  slope  discontinuities  are  less  efficient  in  the 
early  stages  of  loading  a beam,  but  become  more  efficient,  i.e., 
give  closer  values  for  a given  computation  effort,  as  the  yield- 
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point  load  is  approached.  A plate  model  in  which  slope  discontinuities 
provide  the  only  strains  has  been  developed  [37]  and  produced 
some  valid  results.  Work  is  currently  in  progress  [38] 
extending  these  ideas  to  problems  in  plane  strain. 

Even  when  the  rate  problem  has  been  solved,  there  are 
questions  to  be  answered  concerning  the  application  of  Eq.  (11.1). 

A plastic  element  requires  f = 0.  Since  the  yield  surface  f = 0 
is  convex  [23]  , if  the  stress  point  is  not  stationary,  it  must 
move  along  this  convex  surface.  However,  (11.1)  is  motion 
along  a straight  line  tangent  to  the  surface.  Thus  at  time 
tQ  + At,  the  predicted  f will  be  positive  in  violation  of 
(8.6a) . 

A usual  method  of  resolving  this  dilemma  is  to  first  use 
(11.1)  to  predict  a new  a which  will,  in  general,  have  a 
positive  f in  each  plastic  element.  We  then  replace  a by 
Scr  where  the  scalar  S is  chosen  so  that  f(Scr)  = 0 in  each 
plastic  element.  However,  the  stresses  8a  will  not  precisely 
satisfy  the  equilibrium  equations  so  that  some  type  of 
iteration  must  be  performed. 

If  the  yield  function  f is  piecewise  linear,  then  this 
problem  disappears.  The  stress  point  may  lie  on  one  or  more 
linear  faces  so  that  f and  X are  replaced  by  f^  and  X^ 
throughout  (8.6)  and  (8.6e)  is  replaced  by 

^ = E X.  Vf.  (11.2) 

1-1 

If  the  external  loading  is  linear,  Eq.  (11.1)  will  now  be 
exact  until  either  an  elastic  element  becomes  plastic  or  until 
a plastic  element  encounters  a new  face. 

The  yield  function  is,  of  course,  a physical  property  of 
the  material.  However,  it  can  only  be  measured  with  limited 
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accuracy,  so  that  it  is  reasonable  to  approximate  it  by  any 
convenient  function.  For  finite  element  models,  a piecewise 
linear  approximation  has  the  obvious  advantage  of  eliminating 
the  problem  associated  with  (11.1)  and  any  strictly  convex 
yield  function.  A systematic  procedure  for  constructing  a 
piecewise  linear  approximation  to  a given  nonlinear  f is  given 
in  [39]. 

A third  possible  solution  to  this  difficulty  is  to  introduce 
a small  plastic  viscosity  [40,  41].  The  resulting  visco- 
plastic model  permits  expansion  of  the  yield  surface,  but 
for  small  v the  expansion  would  be  severely  limited  so  that 
the  solution  will  be  very  close  to  that  for  a perfectly- 
plastic  material.  Indeed,  for  some  real  materials  the  visco- 
plastic model  may  be  more  realistic  than  the  perf ectly-plastic 
one . 

The  preceding  discussion  of  (11.1)  has  been  implicitly 
based  upon  a model  with  constant  stress  in  each  element.  For 
a more  complex  element,  the  problem  would  be  equally  present, 
and  all  of  the  suggested  remedies  would  be  more  difficult  to 
apply.  For  this  stage  too,  then,  there  are  strong  argximents 
in  favor  of  using  a large  number  of  simple  elements  in  the 
solution  of  plasticity  problems. 

12.  Conclusions . It  is  evident  that  the  preceding  sections 
make  no  claim  to  completeness  either  in  content  or  in  references. 

We  have  tried  merely  to  sample  the  vast  field  of  plastic  structures, 
beginning  with  the  simplest  truss  and  beam  problems  in  the  hope 
of  rousing  the  reader's  interest.  Many  interesting  and 
important  topics  have  had  to  be  omitted  entirely.  Among  these 
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omissions  are  plastic  buckling,  effect  of  geometric  changes, 
and  dyncimic  effects,  and  we  have  not  even  mentioned  the  vast 
literature  on  plastic  shells  [42], 

The  reader  interested  in  more  details  concerning  the 
material  covered  will  find  most  of  it  and  more  in  [3] . For 
the  material  not  covered,  almost  every  new  issue  of  most 
mechanics  journals  contains  articles  extending  the  theory  of 
plasticity  and  the  list  of  problems  which  it  has  solved. 
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(a)  Elastic/strain  hardenine 

(b)  ElasUc 

(c)  Rigid/strain  hardening 

(d)  Elastic/perfectly  plastic 

(e)  Rigid/perfectly  plastic 

Fig.  2— Ideal  stress-strain  curves 
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Fig.  4 — Load-displacement  curves  for  truss  (letters  in 
parentheses  refer  to  the  materials  of  Fig.  2) 
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Fio.  6 . Moment  curvature  relations.  OABD,  perfectly  plastic  material.  OABC, 
mild  steel.  OAED,  ideal  section.  OFED,  rigid  plastic  material. 
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Fio.  8 Collapse  of  simple  frame,  (a)  Ixiaded  frame.  (6)  Collapse  mode  for  beam 
failure  only,  (e)  Collapse  mode  for  p.anel  failure  only,  (d)  Combined  collapse  mode, 
(c)  Moment  distribution  for  (b)  (not  admissible).  (/)  Moment  distribution  for  (e) 
(statically  admissible),  {g)  Moment  distribution  for  (d)  (not  admissible). 
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bic.  7 Collapse  of  singly  indeterminate 

beam,  (a)  Loaded  beam.  (6)  Elastic 
moments,  (c)  One-hinge  moment  dis- 
tribution. (d)  Beam  with  one  hinge, 
(e)  Collapse  moments.  (/)  Collapse 
mechanism. 
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Pig.  23  Stress  and  velocity  profiles  for  annular  plate:  (a)  static- 
ally admissible  stress  profile;  (b)  statically  inadmissible 
stress  profile;  (c)  correct  stress  profile;  (d)  kinematically 
inadmissible  velocity  profile;  (e)  kinematically  admissible 
velocity  profile;  (f)  correct  velocity  profile 


